Fingerprints of intrinsic phase separation: 
magnetically doped two-dimensional electron gas 
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In addition to Anderson and Mott localization, intrinsic phase separation has long been advocated 
as the third fundamental mechanism controlling the doping-driven metal-insulator transitions. In 
electronic system, where charge neutrality precludes global phase separation, it may lead to various 
inhomogeneous states and dramatically affect transport. Here we theoretically predict the precise 
experimental signatures of such phase-separation-driven metal-insulator transitions. We show that 
anomalous transport is expected in an intermediate regime around the transition, displaying very 
strong temperature and magnetic field dependence, but very weak density dependence. Our predic- 
tions find striking agreement with recent experiments on Mn-doped CdTe quantum wells, a system 
where we identify the microscopic origin for intrinsic phase separation. 



Two basic routes to the metal-insulator transition have 
long been known: the Anderson fl'| (disorder-driven) and 
the Mott [2] (interaction-driven) mechanism for electron 
localization. In many novel materials, however, the inter- 
mediate "bad metal" regime is reached by lightly doping 
a parent insulator. Here, a new possibility emerges - 
the tendency of the few available charge carriers to "con- 
dense" into a Fermi liquid - leading to phase separation 
However, charge neutrality precludes global phase 
separation, and the ground state becomes inhomogeneous 
- even in the absence of impurities or disorder. Such 
inhomogeneities can assume various patterns [i5], includ- 
ing bubbles, stripes, or checkerboard states, and may also 
cause glassy behavior even in the limit of very weak dis- 
order [g, 0| ■ The relevance of nano-scale phase separa- 
tion, in its many proposed forms, has been advocated for 
the bad metal regimes of cuprates 0, 01 and manganites 
Q, and even low density two-dimensional electron gases 
(2DEG) 0. Still, in most of these cases, no clear and con- 
clusive evidence has been presented that intrinsic phase 
separation (IPS) - as driving force for the metal-insulator 
crossover - dominates the bad metal regime. While good 
evidence exists that here the carrier concentration is in- 
deed inhomogeneous, this could very well be caused by 
impurity effects. 

To make progress in understanding the relevance of 
IPS, it would be very useful and important to identify a 
model system, that would satisfy the following criteria: 
(1) the microscopic origin for the tendency of phase sepa- 
ration can be established; (2) the system is clean enough 
to avoid extrinsic impurity effects; and (3) transport sig- 
natures can be identified that make it possible to distin- 
guish IPS from spurious disorder effects. In this letter 
we concentrate on a specific physical system - magneti- 
cally doped 2DEG - where all these criteria can be sat- 
isfied. We formulate a simple theoretical model that de- 
scribes IPS in this system, and present theoretical results 
predicting its transport properties in the metal-insulator 
crossover regime. 

Magnetically doped 2DEG: ideal IPS model system. 
The criteria for observing IPS are not easy to satisfy in 
most available systems. In many cases, such as for the 
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Figure 1: (Color online) In presence of magnetically-driven 
IPS, low-density carriers are trapped in ferromagnetic "bub- 
bles" from which they can escape only by thermal activation. 
Such activation-limited transport results in resistivity which 
is a strong function of temperature and magnetic field, but a 
very weak density dependence - the smoking gun for IPS. 

important family of doped cuprates l^], its microscopic 
origin remains highly controversial. In addition, chem- 
ical doping is typically linked with appreciable amount 
of disorder created by dopant ions; such is the situation 
in most manganite oxides and diluted magnetic semicon- 
ductors (DMS). Very recent experimental work fio'l, how- 
ever, has presented striking results on a model system, 
which seems to be an ideal setting for observing IPS. 
Here, a standard 2DEG (formed in CdTe quantum wells) 
has been doped by Manganese - a situation that can be 
described as a two-dimensional DMS system. Because 
this 2DEG is created by (iodine) modulation- doping, the 
impurity ions are far from the mobile electrons, which 
dramatically reduces the disorder effects as compared to 
the well-known CdMnTe bulk (three dimensional) ma- 
terials. As a result, the metallic behavior here persists 
down to remarkably low carrier densities, even under sig- 
nificant Mn doping. In addition, in this materials (as 
opposed to the more familiar GaMnAs), Mn is neither a 
donor nor an acceptor. 

In this way, devices have been fabricated where, around 
the critical carrier density n^, there are only a few per- 
cents of carriers per Mn ion. The possibility of reach- 
ing this regime is very significant, because it is precisely 
where one expects the Mw-induced magnetism to favor 
phase separation [§1 . Indeed, both experimental and the- 
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oretical considerations have established fio'l that CdM- 
nTe is, under sufHcient carrier doping , a typical Zener 
(double-exchange) ferromagnet (ill . . In such materi- 
als sufficient doping induces ferromagnetic (FM) correla- 
tions between the Mn spins, which can prevail over the 
usual Mn-Mn antiferromagnetic (AFM) superexchange 
that exists even in absence of carriers. At very low dop- 
ing levels, however, the situation is more subtle. Here 
the average kinetic energy gain obtained by FM order- 
ing is not sufficient to balance the energy cost of AFM 
superexchange that opposes it. Global FM is then not 
possible, but a compromise is reached through phase sep- 
aration. FM ordering then prevails only in carrier-rich 
regions, where phase separation makes it possible to lo- 
cally achieve a sufficient carrier concentration favoring 
Zener ferromagnetism. 

Model. While the 2D CdMnTe system offers an es- 
pecially attractive venue to experimentally observe it, 
magnetically-induced IPS in a clean low density electron 
gas is a very robust phenomenon [81 , not confined to two- 
dimensional systems or the specific material properties of 
CdMnTe. To emphasize its universal properties we con- 
centrate on the simplest model that illustrates qualitative 
aspects of such an IPS, and is given by the Hamiltonian 

H = H,nag + Hcoul + Himp. (1) 

The ffist term here describes a (random) lattice of classi- 
cal spins Si (representing the Manganese local moments) , 
interacting with conduction electrons through "s-d" ex- 
change interactions Jsdi and with each other through 
(AFM) superexchange interactions J^x- Following An- 
derson and Hasegawa jlH, for simplicity we focus on the 
Jsd ~^ oo limit, where the conduction electron spins com- 
pletely align with the local moments 

H„rag = - ^teos {9ij/2) {c\cj + c]c,) + Jex ^ S.^Sj. 

ij <ij> 

(2) 

Here, 6ij — arccos(S'iS'j) represents the relative angle be- 
tween the neighboring Manganese spins, and cl is the 
conduction electron creation operator corresponding to 
lattice site i. The second term Hcoui represents the 
long-range Coulomb repulsion between conduction elec- 
trons; its principal role is to prevent global phase separa- 
tion in order to keep the charge neutrality of the system. 
The last term Himp describes (weak) disorder scattering, 
which in our system takes the form very similar as in 
other 2DEG materials not containing manganese spins. 

IPS at T=0. To reveal the physical mechanism for 
IPS in our model, and understand the respective roles 
of the two magnetic interactions, we start our analy- 
sis by first constructing the ground state phase diagram 
for H„iag- In absence of superexchange {Jex = 0), the 
ground state for any doping corresponds (l2| to com- 
plete ferromagnetic (FM) polarization of all Mn spins, 
as this allows for maximum kinetic energy of conduction 
electrons. For J^x > 0, however, at infinitesimal carrier 
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Figure 2: (Color online) DMFT phase diagram. Shaded areas 
show the IPS regions at three different values of Jex. As Jex 
increases, the critical density ric increases, while the critical- 
end-point temperature Tc decreases, until IPS is suppressed 
at J > Jex ~ 0.025D; the inset shows the analytical (dashed 
line) low-temperature approximation (see text) for Tc{Jex), as 
compared to the numerical solution (dots). Note the "Pomer- 
anchuk effect" - the re-entrance of the coexistence region, due 
to higher spin entropy of the carrier-poor PM phase. 

doping the corresponding energy gain is not sufficient to 
overcome the cost of superexchange energy opposing spin 
alignment. No ferromagnetism is then possible, instead, 
the system undergoes phase separation between a FM 
carrier-rich phase, and a paramagnetic (PM) flf] carrier- 
poor (carrier-free at T = 0) phase. For our classical-spin 
model, the condition for IPS (phase coexistence) is ob- 
tained by equating the (grand) thermodynamic poten- 
tials of the respective phases, viz ripM = ^pm ■ In the 
FM phase 

p{w)wdw + Jex/'^ ^ ^ / p{w)dw, (3) 

-D J -D 

where p(w) is the density of states of (spin polarized) 
conduction electrons, and D ~ t is the corresponding 
half-bandwidth, while Qpm = [18]. We find that IPS 
emerges when the carrier density n (per AIn ion) satisfies 
the condition n < ric, with 

n,^ji/(^+^), (4) 

and fj,c + D ^ jtx '^^'^^- In the entire IPS regime, both the 
local carrier concentration in the FM phase npM — Uc, 
and the corresponding chemical potential fipM^Pc re- 
main pinned at their critical value (independent of the 
average density n). As a result, the activation energy 
gap to the (PM) conduction band fio'l Eg = Ec — fi also 
remains n-independent. This simple mechanism, com- 
mon to any example of IPS, is the essential cause for 
very weak density dependence of (activated) transport in 
this regime. 

Finite temperature phase diagram - DMFT solution. 
At finite temperatures, the situation is more compli- 
cated. First, a fraction of electrons are transferred from a 
carrier-rich FM domains into into the carrier-poor PM re- 
gions. Second, at finite temperatures the entropy of spin 
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fluctuations strongly favors the PM phase. To quantita- 
tive account for both of these effects, we now construct a 
complete phase diagram for our system using dynamical 
mean- field theory (DMFT). This approach, which is for- 
mally exact in the limit of large coordination, has recently 
been shown Q to provide quantitatively accurate results 
even in low dimensions, for the class of models we con- 
sider. To determine the domain of magnetically-driven 
IPS, we first focus on effects described by Hmag', we will 
later return to the role of Hcoui a-nd its consequences on 
transport. 

Within DMFT, the lattice model is reduced to solv- 
ing a single-site quantum impurity model supplemented 
by an appropriate self-consistency condition. For our 
Anderson-Hasegawa model, the local effective action 
takes the form 

Seffii.S,) = / dT dT'4iT)[SiT~T'){^r+^i) 

Jq Jq 

+ A,{T,T',S)]c,{T')-l3J,,mS„ (5) 

Here, we have used functional integration over Grass- 
mann fields Ci(r) describing electron at site i. The form 
of the self-consistency condition defining the "cavity field" 
Ai(T, r',5) depends on the electronic dispersion [l3| . 
Within the DMFT, however, the results do not quali- 
tatively depend on the precise form of the band struc- 
ture or the dimensionality; without loss of generality we 
therefore use a simple semi-circular model bare density 
of states p(){uj) — j^y/D'^ — uP-. Using standard methods 
[III, we find 



-G«H-5), (6) 
and (Ml{wn) 



where Ga2(w„) =< c'l'(w„)c(w„) > 
Sc^ (wn)c{^w.m) >Scff ^^c the local Green's functions, and 
m —< S >s^ff is the magnetization. For simplicity, we 
also use Ising {S = ±1), as opposed to Heisenberg spins; 
we checked that within DMFT this only leads to very 
small quantitative changes. These DMFT equations are 
easy to solve numerically as a function of the carrier den- 
sity and temperature, for several values of the superex- 
change J^x'-, the resulting phase diagram is presented in 
Fig. 2. 

Temperature- driven IPS - Pomeranchuk effect. Our 
numerical solution reproduces Eq. [3]at T = 0, but inter- 
esting re-entrance (Fig. 2) is found at T > 0. Physically, 
the temperature increase favors phase separation, be- 
cause of the larger spin entropy Sspin ~ kBln2 associated 
with the carrier-poor PM phase, very similar to the ef- 
fect Pomeranchuk predicted for He^ solidification. Very 
similar Pomeranchuk-like behavior of phase separation 
behavior has been discussed for other systems Q, and 
generally it may give rise to pronounced resistivity max- 
ima just on the metallic side of the MIT. To clarify this 
re-entrance, we note that the leading low-temperature 
behavior can be analytically obtained from Eqs. (j5l6p . 
by performing an appropriate Sommerfcld expansion. To 
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= 0) -h 0(T2), while 
giving an expression 



for the high-density boundary of the coexistence dome 
(dashed line on Fig. 2) of the form 
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(7) 



A similar argument gives an expression for critical tem- 
perature of IPS 



(8) 



where for our model J^^ « 0.0251), and C « 0.88. This 
shows, in perfect agreement with our numerics (see inset 
of Fig. 2), that while ric increases with J^x (see Eq. U]), 
the critical temperature Tc decreases, until phase sepa- 
ration is suppressed for J^x > Jex- 
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Figure 3: (Color online) (a) and (b) show the resistivity as 
function of temperature at different carrier densities. For 
comparison, we contrast the experimental results [To| | (left 
panel) and our DMFT predictions (right panel). A sharp up- 
turn in resistivity in the low temperature regime, and its weak 
density dependence are direct consequences of IPS. Plots (c) 
and (d) show the same data as function of carrier density at 
fixed temperatures to emphasize weak density dependence of 
resistivity below the percolation threshold. It also reflects the 
fact that within this region activation energy as well as the 
size of the PS regions is density independent (see text). 

Phase separation and transport. As density increases 
in the IPS region, the volume fraction x of the carrier-rich 
EM phase determined from n = xnpM + (1 ~ x)npM also 
increases, until a percolation concentration is reached 
and metallic transport sets in. In presence of long- 
range Coulomb interactions, the characteristic size of 
EM) carrier-rich domains is limited by the charging en- 
ergy, which can be easily estimated. 

However, within a commonly-used random-resistor pic- 
ture of transport (which is well known to be adequate for 
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the kind of system we consider) the critical behavior is 
largely insensitive to the individual domain size - with 
the only relevant parameter being the relative volume 
fraction x. Thus, to obtain a quantitative estimate of 
transport properties for our model, it sufHces to calcu- 
late the local temperature-dependent resistivity for each 
of the two phases PFAf/PA/ ('^i T). This task is easily car- 
ried out using standard methods [13, 14] within DMFT 
model allowing us to determine the resistivity from the 
solution of our DMFT equations. This calculation ac- 
counts for density and temperature dependent resistivity 
PFM/PMin,T) due to scattering of conduction electrons 
off Manganese spins in the respective phases. 

To make a meaningful comparison with experiments, 
however, we also have to account for conventional im- 
purity scattering while considering the total resistivity 
behavior. This weak impurity scattering in our system 
which, by the way, is unrelated to the physics of IPS 
and is more pronounced at higher temperatures, is of 
the same form as in conventional 2DEG systems [T5| . 
In particular, as can be seen from the high-temperature 
transport data of Ref. [IQj, it has very weak tempera- 
ture dependence but an appreciable density dependence, 
which we obtain by fitting the experimental data, giving 
Pimp ~ Po exp{-A?i} + C, (here we used po = 3600pMott, 
A = 850 and C = Q.2pMotu resistivity is in units of pMott 

The net resistivity within the IPS region is then calcu- 
lated by solving the random resistor problem, where we 
use the familiar effective-medium approximation (with 
Pmed being determined as in (Eq. 5.7) of |16]), giving 

p{n,T) = Pmed{Pfm,Pp7n) + Pimp- 

The resulting resistivity curves are given in Fig. 3; as 
we can see, our model seems to capture all the qual- 
itative features of the experimental data. Below the 
(temperature-dependent) percolation threshold, where 
we see a sharp increase in resistivity, the transport is 
activated with essentially density-independent activation 
energy. At higher densities, we note the well-pronounced 
resistivity maximum at low temperatures, reflecting the 
re-entrant Pomeranchuk-like IPS behavior we described 
above. 

CMR behavior. As in other double-exchange ferromag- 
nets, a moderate magnetic field reduces the IPS region, 
quickly lowers the relevant activation energy, and leads 
to colossal magnet-resistance (CMR) behavior (Fig. 4). 




Carrier Density Temperature 

Figure 4: (Color online) Phase diagram (a) and transport (b) 
in presence of an external magnetic field h. Panel (a) shows 
how the temperature and density parameter range of IPS is 
quickly reduced in nonzero magnetic field, and (b) depicts the 
negative magnetoresistance within this IPS region (also seen 
in [T9|). The inset of (b) shows the decrease of activation 
energy Eg (which was extrapolated from the resistivity vs. 
temperature plots) with magnetic field. 

However, this behavior is restricted to the coexistence 
region, and as such can serve as an indication for IPS 
behavior-again in close agreement with experiments[10]. 

Conclusions. In this letter we presented and solved 
a simple microscopic model for magnetically-driven IPS 
behavior. Given the simplicity of our theoretical model, 
the agreement with all the qualitative features of the ex- 
perimental data is remarkable. Most importantly, we 
observe how, within the IPS regime, all the resistivity 
curves "merge" together below a temperature scale set 
by Tc for phase separation. This striking feature simply 
reflects the fact that the chemical potential - setting the 
activation scale in this problem - is fixed and density- 
independent everywhere under the IPS dome. Such be- 
havior is quantitatively reproduced in our model for mag- 
netically driven IPS. However, very general thermody- 
namic principles guarantee that the very same behavior 
must be at play in virtually any system where transport 
is driven by a tendency to inherent phase separation - 
and not the disorder effects. In this work we identified 
the first physical system which presents plausible and 
convincing evidence for the IPS scenario. From a more 
general point of view, our results indicate how should one 
search for other examples where this fascinating physics 
may be recognized and properly interpreted. 
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